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TECHNICAL NOTE 5357 


THE EFFECTS OF VARIOUS PARAMETERS, INCLUDING MACH 
NUMBER, ON PROPELLER-BLADE FLUTTER WITH 
EMPHASIS ON STALL FLUTTER^ 

By John E. Baker 

SUMMARY 


The effects of many of the paxameters significant to vring flutter 
were studied experimentally on several untwisted rotating models to 
determine their significance with respect to stall flutter of propeller 
"blades. The parameters included torsional stiffness, section thic kn ess 
ratio, sweepback, length- chord ratio, section center- of-gravity location, 
"blade taper, Mach number, and fluid densiiy. In order to check on the 
effects of blade twist, one model which had spanwise twist was studied. 

The blade angles of the models were generally varied ftom low values to 
beyond the stall. 

The experimental results for the flutter speed are presented in the 
form of flutter-speed coefficients ('7/bc‘tt,)o.8L? "where the quantities V 
and b are the resultant velocity and semichord, respectively, taken at 
0.8 blade length L, and c%, is the natural first- torsion circular fre- 
quency of the blade. The minimum values of this flutter-speed coefficient 
were found to be slightly greater than 1.0 for subcritlcal Mach numbers. 
The parameters that produced a significant increase of this flutter-speed 
coefficient were forward movement of the section center-of-gravlty 
location, sweepback, increase of the section thickness ratio^ and Mach 
number at supercritical speeds . In order to maintain satisfactory aero- 
dynamic efficiency at high speeds, however, thin blade sections are 
required. The largest effects on the flutter- speed coefficient were 
produced by the section center-of-gravity location and the Mach number. 

The effect of Mach nimiber was of such significance that it is possible 
to present a tentative cariterlon for designing completely flutter-free 
thin propeller blades to operate at supersonic and supercarltical speeds . 
This criterion is given by the design parameter (13 CH[,/c)q^qq (where the 
quantity c is the sound speed of the operating medium) and the present 
investigation indicates that propeller blades having values of this 
parameter above 0.50 should be entirely free of flutter. It is also 

^-Supersedes the recently declassified NACA RM L50L12b, "The Effects 
of Various Parameters Including Mach Nimiber on Propeller-Blade Flutter 
With Emphasis on Stall Flutter" "by John E. Baker, 1950. 
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possible to consider various operating procediires for propellers not 
satisfying the criterion as a means for avoiding flutter. 


INTRODUCTION 


The propeller- flutter problem is an old one. Papers have been 
written on the subject from World War I to. the present, n.nfl some of the 
significant results are summarized in reference 1. In the past, propeller 
flutter has been studied primarily by trial-and-error methods on specific 
propellers. Therefore, much confusion and contradictory data exist at 
the present time. 

In general, flutter can be avoided by making the blade sections 
sufficiently thick; however, recent aerodynamic investigations have shown 
that propellers for futvire high-speed aircraft must be thin in order to 
obtain satisfactory performance. Thus, propeller flutter, which has been 
considered chiefly as a nuisance In the past, now becomes a critical 
design problem. The trial-and-error methods of the past have become 
Inadequate for the design of thin supersonic and transonic propellers. 
Accordingly, a fairly comprehensive experlj^ntal study has been made to 
determine the effects of various parameters on propeller flutter. The 
parameters studied include torsional stiffness, section thickness ratio, 
sweepback, length-chord ratio, section center -of -gravity location, blade 
taper, blade twist, Mach number, and density of the operating medium. 

Blade angles were generally varied from low lift to beyond the stall. 

Propeller flutter, as described, for example, in references 2 and 3 > 
can be separated into two main types, clas$ical flutter and stall flutter. 
Classical flutter exists at low angles of attack where the flow can be 
considered potential and, hence, the aerodynamic forces can be evaluated 
by the use of potential-flow theory. Stall flutter is encoimtered at 
high angles of attack where the flow is brqken down and potential theory 
falls to apply. Classical flutter is generally characterized by a coupling 
of the bending and torsion vibration modes; whereas stall flutter occurs 
primarily in the torsion mode. The classical-flutter frequency usually 
falls between the first-torsion and first-bending frequencies, but the 
stall-flutter frequency is nearly the same as the natural first-torsion 
frequency of the blade. There is no definite boundary between stall 
flutter and classical flutter, and a continuous merging exists. The 
natural phenomena Involved in this merging : are st ill uncertain although 
various attempts have been made to associate them with the behavior of 
the static lift curve. Stall-flutter speeds have been found to be lower 
than classical-flutter speeds. 

The designer is primarily Interested in being able to calculate the 
propeller-flutter conditions in connection Iwith the possible operating 
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conditions. At the present time, no theories are established that can 
predict adequately stall-flutter speeds, hut scone theories exist that 
apply to classical flutter of propellers, (references 4 to 6) and a 
brief discussion of these theories is given. 

Propellers under normal flight conditions generally operate with 
the blade sections at low angles of attack and would be subject to clas- 
sical flutter. Under these operating conditions, the flutter speed is 
high, and generally there is an appreciable margin of safety between the 
operating speed and the flutter speed; however, during the take-off 
period, the propeller blade sections may operate at high angles of attack 
and would be subject to stall flutter. Under these operating conditions, 
the flutter speed is generally quite low, and the propellers must operate 
near the flutter condition where the margin of safety is fairly small. 

As a consequence, propellers are required to operate in the stalled con- 
dition without dangerous flutter. Since stall flutter is the critical 
design condition, and because no adequate theories are established for 
predicting stall-flutter speeds for propellers, the investigation was 
devoted primarily to obtaining information on stall flutter, and, in 
particular, the mlniTnuTn stall-flutter condition. 

SYMBOLS 


potential nonsteady aerodynamic coefficients 


nondimens ional elastic axis position referred to semichord 
measured from midchord 

blade semichord, feet 

blade semichord at reference station, feet 
sound speed of operating medium, feet per second 
section lift coefficient 

section center-of -gravity location, percent chord 
shear modulus of elasticity, pounds per foot^ 
torsional stiffness, pound-feet^ 
elastic axis location, percent chord 


•^ch 

Aca 

Aah 

■^aa 




b 

br 

c 

Cl 

CG 

G 

GJ 

EA 
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"blade first-torsion frequency*, cycles per second 

"blade first -"bending frequency, cycles per second 

■blade-"bendlng-deflection function in terms of tip deflection 

"blade-torslonal-deflection fuhction in terms of tip deflection 

structural danqping coefficient in bending as used in 
reference 7 

structural damping coefficient in torsion as used in 
reference 7 

bending deflection of blade, feet 
hub radius, feet 

polar* moment of inertia about elastic axis per unit length, 
slug-feet^ per foot 

polar moment of inertia about the section center-of-gravity 
location per \init length, slug-feet^ per foot 

reduced frequency ("bcu/'V) 

blade length, feet 

Mach number 

blade mass per unit length, s3,ugE per foot 

nondimensional radius of gyration of blade section about 
elastic axis (l^ymb^) 

nondimensional radius of gyration of blade section about the 
section center-of- gravity location (l 0 Q/m"b 2 ) 

section thickness, feet 

resultant velocity, feet per second 

distance from blade root, feet*. 

nondimensional center-of -gravity position measured from 
elastic axis in terms of semichord 


blade angle, degrees 



MCA. TN 3357 


5 


n* 

6 

K 

p 

G 

0 


mass density of blade material, slugs per cubic foot 
nondimens ional distance from blade root (x/L) 
nondimensional distance from center of rotation 


X + H 
H + L 


torsional deflection of blade, radians 
blade mass-density ratio (jtpb2/m) 

mass density of operating medium, slugs per cubic foot 
blade solidity at 0.8 l (2b/2jr(0-8L + H) ) 

aerodynamic helix angle, degrees 


CD 

blade 

circular 


blade 

circular 

per 

second 


blade 

circular 

per 

second 


APPAEATUS AND TEST METHODS 


The apparatus and testing techniques herein described are similar 
to those used in the tests of reference 2. The propeller models were 
operated in the Langley vacuum sphere in which provisions are made for 
operating in air or Freon- 12. Freon-12 is a convenient medium for 
studying the effects of Mach number because its sound speed is about 
500 feet per second at room temperature (reference 8) . The propeller 
models were rotated by means of a 500 -horsepower electric motor (fig. 1) 
and operated at zero forward velocity except for Induced flow. Bending 
and torsion oscillations of the blade were recoixied by an oscillograph 
with the aid of wire strain gages on the blade (see sample record in 
fig. 2) . The rotational speed was also recorded on the same record, 
which, for zero forward velocity. Is equivalent to the resultant velocity. 
A few total-pressure measurements were obtained in the wake by means of 
a survey rake located about 0. 17 propeller diameter behind the propeller 
disk. 


Flutter runs were generally made in air at l/h, 1/2, and 1.0 
atmosphere pressure, but only the data obtained at 1.0 atmosphere are 
presented herein, with the exception of data for the studies of density 
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and Mach number. Blade angles were usually varied from low lift to 
beyond the stalling angle. Flutter was not studied at zero thrust on the 
\mtwlsted blades because wake flutter, such as that described In refer- 
ence 3 , occurred. During each flutter run, the rotational speed was 
gradually Increased until flutter was obse'rved, at which point a record ' 
was taken. A few attempts were made to go thixmgh the flutter region 
at stall, but the flutter was usually too -violent to do so. 

The effect of Mach n\amber was studied by operating In various 
mixtures of air and Freon- 12 In order to Vary the so\md speed of the 
operating medium. This technique made It possible to obtain a range of 
Mach numbers at any given rotational speed. The density was held constant 
at about 0.0011 slug/ cubic foot for the veirlous mixtures by varying the 
pressure of the operating medium. ' 

The flutter models with their Identifying numerical designations 
and their significant parameters are listed In table I. The parameters 
studied, the range of values covered, and the’ models used to study them 
are described In table II. Information about the blades, which Is not 
listed In tables I and II, Is described as follows: 

(1) The swept back models were swept from a radial line with the 
sweepback beginning at the root of the blade, as indicated by the dashed 
outline In figure 1. 

(2) Models la and lb were successively shortened to change the 

length-chord ratio. i: 

(3) The section center-of-gravlty location was varied by the use of 
different blades with brass Inserts cycle-^lded In the blades near the 
leading edge so that . the section contour r^ained unaffected. 

(4) Model 2 was nearly Identical to model la and was twisted 
manually to beyond the yield stress, resulting In a set twist of 17 ° at 
the tip, and the angle of twist varied linharly along the span. 


RESUITS AND DISCUSSION 
Considerations on Method of Presentation 

Reference section .- The experimental’ data are presented showing 
the effects of the various parameters studied on the flutter-speed coef- 
ficients. The data shown are all referred; to the 0. 8-blade- length 
position which, for propellers having largfe hub diameters, would result 
in a more representative reference section: than would result If a 
standard radius location were used. 
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Lift coefficient . - The blade angles shown are the blade-angle set- 
tings at 0. 8 blade length referred to the plane of rotation. The relation 
between blade-angle setting and lift coefficient is distorted because of 
the effects of induced flow and blade twisting due to centrifugal forces 
and aerodynamic forces. Since the designer is primarily concerned with 
lift coefficient, some of the wahe-survey data were evaluated to yield 
lift coefficients. The individual wahe surveys are Incomplete and, as 
a consequence, considerable scatter of the pressure measurements is 
present; however, as a matter of Interest, a band showing the appivDxi- 
mate values of experimental lift coefficients corresponding to the various 
blade-angle settings is shown in figure 3- This figure is applicable 
in general to the models having 0. 333-foot chord and blade length of 
1.788 feet, with which most of the low Mach number data were obtained. 


Flutter-speed coefficient . - The flutter velocity is a function of a 
great number of parameters: 


— = 
bo^ I 


Oiu rs 1, 


The ratio V/boi^ is designated as the flutter-speed coefficient, the 
value of which is dependent on the large number of parameters. This 
coefficient is taken at a reference section which is 0.8 l for the data 
shown herein. The purpose of this investigation is to determine the 
effects of many of these parameters on the flutter-speed coefficient. 

Before discussing the effect of the parameters studied on this coefficient, 
it appears advisable to point out the significance of the flutter-speed 
coefficient and its component parts in order to interpret correctly the 
applicability of the data presented herein. For comparison purposes, 
assume that a certain flutter- speed coefficient is given, in other words. 


t: — = Constant 


For this condition, an increase in the semlchord is accompanied by a 
proportional increase in the flutter speed provided the torsional 
frequency remains constant. 


The semlchord can be varied without changing the torsional frequency 
if the airfoil section is \mchanged, as is illustrated by considering the 
first-torsion frequency equation for a uniform beam: 



J ~ A(2b)t8 where A is a constant 
B(j.t(2b)8 where B is a constant 
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Then 

jt_ / QA(2b)t3 ^ Jt_ M 
2LyB^t(2b)3 2b2LVBn 


( 1 ) 


For a given value of the flutter-speed coefficient, increasing the 
torsion frequency will be accompanied by a proportional rise in the 
flutter velocity if the chord at the reference section is held constant. 

It can be seen in equation (l) that tw^o m^hods of raising the torsional 
frequency are increasing the section thickness ratio and decreasing the 
blade length. Reference 9 indicates that tapering the blade chord will 
also raise the torsional frequency. 

Appropriate care should thus be e;xercised in the interpretation, in 
terns of actual flutter speeds, of results which are presented as flutter- 
speed coefficients. 


Experimental Data end Discussion 

The parameters studied and the figures in which the data are pre- 
sented are listed in table II. In the experimental investigation, the 
various parameters were isolated where: possible, and, correspondingly, 
the data showing the effects of each parameter are presented in separate 
figures. In figures 4 to 11, the oordlhate is the flutter- speed coeffi- 
cient (V/bu^)Q^0L an<3. the abscissa is blade angle Po.8!L param- 

eters studied are torsional stiffness (fig. 4), blade taper (fig. 5)^ 
blade twist (fig. 6), length-chord ratio (fig. 7)^ density of the 
operating me dlinn (fig. 8), section thickness ratio (fig. 9 ), sweepback 
(fig. 10), and section center-of -gravity location (fig. 11). The effects 
of Mach number are shown in figure 12,' where flutter-speed coefficients 
for a given blade angle are plotted as a function of Mach number. These 
parameters are discussed in this section. It is noted that many of the 
flutter curves are not completely filled in at low blade angles. For 
these cases, the flutter speed has become, higher than the maximum safe 
operating speed of the blades. 

The flutter data given in figures; 4 to 11 were obtained under condi- 
tions of subcritlcal flow, that is, with gubcritlcal operating speeds 
at the reference section. A significant observation can be made frraa a 
study of the minimum values of the flutter-speed coefficients that occur 
for each parameter studied; namely, thp lowest value obtained for each 
parameter is slightly greater than 1.0. Deviations from this value are, 
therefcre, used as a basis of comparison for variations of each parameter. 

Parameters having little effect oh the minimum flutter- speed coeffi - 
ients . - The parameters that produced no significant Increase of the 
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minimum flutter-speed coefficients are torsional stiffness, blade taper, 
blade twist, length-chord ratio, and density of the operating medium 
(figs, k- to 8). As has been pointed out, the minimum flutter-speed 
coefficients may be unaffected by changing a given parameter, but the 
product bq^ should be examined to determine the effect of the changes 

on the flutter speed. 

In contrast to the insignificant effect of torsional stiffness on 
the TnlniTnuTTi flutter-speed coefficient, a large effect on the flutter-speed 
coefficient at low blade angles is indicated by the data in figure if. 

This effect is in accord with the theory for classical flutter. 

Section thickness ratio .- Increasing the section thickness ratio 
is shown to have seme effect on the Tn-tniminn flutter-speed coefficients by 
the data in figure 9. Increase of the section thickness ratio from 6 to 
9 percent chord raised the minimum flutter-speed coefficient about 
20 percent; however, thick blade sections are associated with greater 
reductions in aerodynamic efficiency at transonic speeds. 

Sweepback . - The flutter data in figure 10 indicate that moderate 
amounts of sweepback raised the minimum flutter-speed coefficient about 
30 percent. In view of the serious structural problems associated with 
sweptback propeller blades, this moderate rise in minimum flutter-speed 
coefficient does not appear to be of 3oiuch practical significance. 

Section center-of -gravity location .- A pronounced effect of section 
center-of-gravlty location is indicated by the flutter data in figure 11. 
Forward movement of the section center-of-gravlty from 48.5 to 37* 4- percent 
chord resulted in a rise of the tti 1 TilTmrm flutter-speed coefficient of about 
60 percent. At 34-. 0 percent chord, the minimum flutter-speed coefficient 
was about 80 percent higher than that for the section center-of-gravity at 
48. 5 percent chord. 

This favorable effect of forward movement of the section center-of- 
gravity location cannot be utilized to a great extent for solid blades 
but, for built-up or hollow sections, some forward movement of the section 
center-of-gravlty location can be realized. However, forward movement 
of the section center-of-gravlty location for operation at supersonic 
speeds may result in some unfavorable conditions. For example, centri- 
fugal force causes the effective elastic axis of propeller blades to 
approach the section center-of-gravlty location. The aerodynamic center 
of pressure is shifted from the subcritical value of quarter chord to 
about midchord at supersonic speeds. If the section center of gravity 
is located far forward, the aerodynamic pitching moment about the section 
center-of-gravity location at supersonic speeds would become negative. 

This negative pitching moment would then add to, rather than oppose, the 
negative pitching moment due to centrifugal force, probably resulting in 
excessive torsional deflections. 
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The data In figure 11 Indicate that forward movement of the section 
center-of-gravlty location has an extremely gi^at effect on the flutter- 
speed coefficients at low blade angles. T^s effect Is to be expected 
from classical- flutter theory. 

Mach number . - The effect of Mach number on the minimum flutter- 
speed coefficients Is beneficial, as Is Indicated by the data In 
figure 12 for two blades each at a constant blade angle. The blade 
angles were chosen to be the angles at which the minimum flutter-speed 
coefficients were obtained on each blade, as shown In figure 5- In 
figure 12, the coefficients remain nearly constant at about 1.1 to 1.2 
up to the vicinity of the critical Mach number at the reference section. 
Further Increases of Mach number result In a rapid rise of the flutter- 
speed coefficients. 

It Is to be noted In figure 12(b) that flutter was encountered at 
several points In the supposedly stable region at a (v/b<j^)o. 8 l 1*7 
at a rotational frequency of one-eighth thp blade torsional frequency. 

The oscillation encountered Is very likely’ caused by strut Interference 
since there are four struts supporting the ’motor. Further Indication 
of Interference Is supplied by the fact that the range of speeds at 
which these oscillations were observed Is Very narrow. 

The significance of the Influence of Mach number Is better Illus- 
trated by replottlng the experimental flutter curve In figure 12(b) In 
the form shown In figure 13. If^ In the f:|.utter-speed coefficient, 
both the ntmerator and denominator are divided by the speed of sound, 
the two nondlmenslonal coefficients, Mach ipimber and bo^c are obtained, 
both taken at 0.8 blade length. These quantities are used as ordinate 
and abscissa In figure 13. Straight lines .radiating from the origin 
indicate constant flutter-speed coefficients. The value of (bo:j^c)Q qj ^ 

at which the turning point of the flutter curve occurs Is considered 
to be of fairly general significance. This conclusion is confirmed 
by test points obtained frcmi whirl tests of full-scale propellers made 
at Wrlght-Patterson Air Force Base and by different manufacturers. 

Since the experimental flutter curve in flgiure 13 is for the blade angle 
at which the mlnlm\3iti flutter-speed coefficient occurred, data at either 
lower or higher blade angles should fall above and to the left of the 
given Instability curve. The portion of the instability curve above the 
turning point could not be investigated with the apparatus available 
for these studies since the flutter encountered was too severe. A given 
propeller would operate on a vertical line designated by a constant value 
of (hah/c)o,8L ^ fixed speed of sound. It can be seen that, for 

blades having low values of and operating at the stall 

condition, this line would intersect the flutter curve before supersonic 
speeds are reached, and the blades would experience flutter. However, 
it may be possible to design satisfactory thin propeller blades 
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with enough to permit operation into the supersonic 

speed range without intersecting the flutter boundary. 


Possible Applications 

A design criterion .- A tentative design criterion based on these 
data can be determined and indicates that propellers having values of 
the design parameter (bt%/c)o.8L greater than 0,50 should be entirely 
free of flutter. Many current propellers giving satisfactory service at 
tip Mach numbers near 1.0 have values of the design parameter near 0.40, 
These propellers may flutter at the stall, but whirl tests established 
any flutter which may have been encountered as nondestructive. The value 
of 0.50 is used for the criterion presented because thin blades probably 
could not endure flutter without the danger of fatigue. 

Some blade configurations based on the given design criterion are 
shown in figure lif-. Two designs of constant thickness ratio are shown, 
although structurally this condition may not be too practical. Another 
blade having taper in thickness ratio and constant chord, which may be 
more acceptable, is aJ.so shown. These blade configurations may not be 
ideal in some respects, but it appears possible to construct supersonic- 
type propellers with (b<%/c)o.8L greater than O.50 and, consequently, 
to be completely free of flutter. 

A cycling process .- Many of the super sonic -type experimental 
propellers being considered at the present time have values of (bo:t(/c)o.8L 
of the order of 0.10 to 0.20. It can be seen in figure 13 that such 
propellers would flutter if attempts were made to accelerate them to 
supersonic speeds at the stall condition. There is a possibility that 
these propellers can still be operated at supersonic speeds at stall 
without flutter if they are brought up to speed in a manner to be described. 

The f lut ter -r speed coefficient at seme blade angle lower than the 
stalling blade angle would be greater and would appear in figure 13 as 
a line from the origin of greater slope. The lower lift coefficient would 
raise the critical Mach n'umber, and thus the flutter curve at some 
unstalled blade angle should be similar to the dashed curve in figure 13 . 

The experimental flutter curve is extended in the direction it might be 
expected to go by the dotted 3J.ne. A propeller having (b<xtt/c)o.8L of 

say 0.1^ would intersect with the flutter curve if attempts were made to 
bring it up to supersonic speeds at 2CP blade angle; however, it could 
be accelerated to supersonic speeds at the lower blade angle without 
fluttering. Once the propeller is up to speed, the blade angle could 
be increased to 20^ without experiencing flutter since this condition 
would be above the upper limit of the flutter boundary. It is necessary, 
however, that the operating speed is not reduced enough to Intersect with 
the flutter curve due to the Increased power loading. The reverse of this 
operating cycle would have to be followed in stopping the propeller if 
flutter is to be avoided. 
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This so-called cycling process necessitjhes close attention to 
other parameters which are critical to the low angle-of -attack classical- 
flutter speeds, and involves primarily tohsional stiffness. The data 
in figure 4 demonstrate the great effect of torsional stiffness on the 
flutter speed at low hlade angles. A critical condition for successful 
cycling is that the classical- flutter spejed s appreciably higher than 
the maximum operating speed. This condlt,ion exists when (tx%/c^0.8L 

for a given propeller does not Intersect the flutter curve for the blade 
angles at which the propeller is brought .up io speed. 

Operation above the flutter boundary .- The flutter at the min im um of 
the low-speed flutter curves was generall|y much less violent than the 
flutter at lower blade angles. Attempts were made to operate some of the 
test models into the flutter region. At blade angles corresponding to 
the minimum of the flutter curve, some of| the blades were operated 
successfully without dangerous flutter at higher speeds than those 
indicated by the flutter curvej however, 'the density of the operating 
medium usually had to be reduced considerably before successful operation 
resulted. At lower blade angles, the flatter region could not be 
penetrated without the flutter becoming very severe. 


Comparison of Experiment with GjLasjlcal-Flutter Theory 

As previously discussed, propeller flutter can be separated into 
two main types, classical flutter which occurs at low angles of attack 
and stall flutter which is associated with hl^ angles of attack and 
which occurs at lower speeds than does classical flutter. Since the 
designer is Interested in being able to predict flutter speeds, a survey 
of existing theoretical techniques is desirable. 

At present, no theories are established_that can adequately predict 
stall-flutter speeds for propellers. However, in order to make effective 
use of cycling procedures, knowledge of classical- flutter speeds is 
desirable, so some of the available classic^ propeller-flutter theories 
( references 4 to 6) will be discussed briefly. The theory of reference 4 
uses a differential-equation approach shnllar to that used in wing- flutter 
theory, but, in addition, introduces centrifugal force and moment into the 
equations. Reference 5 uses the same attack to the problem, but with 
more simplifying assumptions which eases ' m;imerical application somewhat. 

The theory of reference 6 utilizes known wing-flutter theory in a manner 
similar to references 4 and 5- The effebt of centrifugal force is included 
in the bending mode, but neglected in the torsion mode. Classical two- 
dimensional oscillating air forces are ubed in all three theories, and 
reference 6 has provisions for using either, campresslble or incompressible 
values. Some computations have been madb in order to ccmpare theoretical 
with experimental results presented herein. The theories referred to are 
quite difficult to adapt to numerical calculj.tlons and generally require “ 
considerable computing time. The theory! of reference however, with 
certain modifications, was used to compute one case. 
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The theory of reference 5 vas developed for application to helicopter 
rotors with the assumption that the root of the "blade is located at the 
center of rotation. This assumption does not lead to great errors when 
applied to helicopters because the huh diameter is generally small with 
respect to the rotor diameter. Since propellers liave much larger hubs, 
the theory of reference 5 bad to be modified to make use of a hub radius, 
which may be as much as 30 percent of the propeller radius. The modifl^ 
theory was used to compute the classical-flutter speed of model k, and 
the result is shown in column (l) of table III. 

Since the existing propeller-flutter theories are quite cumbersome, 
a classical wing- flutter theory (reference 10) was modified to apply to 
propellers. This modification was accomplished by allowing the aerodynamic 
forces to vary along the blade and applying centrifugal-force corrections 
to the static fir st -bending and first-torsion frequencies. This method 
of analysis is discussed in the appendix. The dynamic deflection curves 
were assumed to be the same as for the static case. This method was 
used to compute classical-flutter speeds for three of the models used in 
the current tests, and the results are shown in table HI column (2). 

A comparison of theoretical values in column (2) of table HI with 
experimental results in columns (5) and (6) of table HI shows that 
theoretical predictions are sll^tly lower than the experimental classical- 
flutter speeds, but are possibly adequate for predicting classical flutter; 
however, the theoretical values are considerably higher than the experi- 
mental stall-flutter speeds, which indicate that classical theory, using 
oscillating air forces derived from potential flow, is wholly Inadequate 
for predicting stall-flutter speeds. 

It would be less time consuming to ccmpute the classical- flutter 
speed of a given propeller if two-dimensional wing-flutter theory, rather 
than the more tedious propeHer-flutter theory, could be used. This could 
be done if a representative section on the propeHer blade were established 
at which a flutter speed computed "by two-dimensional theory could be 
applied. Calculations "by the wing-flutter theory of reference 7 were made 
on the three models used to compare theozy and experiment, and the results 
are shown in column (3)^ table HI. On the basis of comparing the two- 
dimensional calculations with the propeller calculations in colxomn (2), 
theoretically derived representative sections are determined and are 
listed in column (4). These results show that a value of 75 percent blade 
length may be adequate for the representative section. 
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CONCLUSIONS 


The effect of many parameters significant to wing flutter as well 
as blade twist was studied on several urltwisted rotating models to deter- 
mine their significance with respect to ;prdpeller stall flutter . The 
experimental propeller- flutter data indicate the following conclusions': 

1. The minimum flutter-speed coefficients obtained at low Mach 
numbers were slightly greater than 1.0. ' 

2 . Forward movement of the section center-of-gravity location, 
increasing section thickness ratio, sweepbaCk, and Mach number at super- 
critical speeds were the only parameters studied that raised the mini- 
mum flutter- speed coefficients appreciably above 1.0. Section center- 
of-gravity location and Mach number appeared to show the most significant 
Increases . 

3. The beneficial effect of Mach n\mber indicates a design param- 

eter which is designated by (^3o^/c)Q^gg|^, where b is the semichord 
taken at 0.8 blade length L, is the natural first-torsion circu- 

lar frequency of the blade, and c is the .sound speed, of the operating 
medium. It appears that a tentative design criterion can be given which 
states that propeller blad.es having (hhfcc./^0.8L than O. 5 O 

should be entirely free of flutter. . 

4. Practical supersonic propellers ; having thin blade sections may 
not satisfy the criterion. A proper cycling procedure would then 
probably be necessary whereby the propeller could be accelerated to 
supersonic speeds at low blade angles. 'For this procedure to be success- 
ful, the class leal- flutter speed must be appreciably hl^er than the 
desired operating speed. Once the propeller is up to speed, the blade 
angle can be Increased to the desired loading conditions without 
encountering flutter. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 12, 1950- 
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APPENDIX 

METHOD OF ANALYSIS 


The classical wing-flutter theory of reference 10 was modified for 
application to propellers in the following manner. 

The equations of equilibrium in the torsional and bending degrees 
of freedom are written in reference 10 with the sweepback terms neglected 
as 


where 


A2 


^M2 + rtpbr^o)^ = 0 

(1) 

^^2 + £E2^ icpbj.^0:^ = 0 

(2) 

&)% • “Jj 



Pl.O/ v2 r -I 2 

4 ^ch[Fh(Tl)J dq 

(3a) 


° ■ Aca)[Fn(il)]l?e(>l)]i>I (3i) 

”2 - [FeC^OAir, (3c) 



i6 


MCA TN 5557 



The air forces are given hy 

Aah - » + (i + * l(i * “)(- f ) ° ■ 2 -(I 

Aaa - - I - - (I + - '‘^)f * 


(3d) 


The quantities F and G are the, real and im^inary parts, 
respectively, of the complex function C = C(k) = F(k) + which is 

associated with the wake and was developed by Theodorsen in reference 11, 

The border-line condition of flutter separating the damped and 
undamped oscillations is determined by a nontrivial solution of the 
homogeneous equations (l) and (2). The, flutter condition is solved by 
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means of the vanishing determinant of the coefficients of the bending 
and torsional motions. 


A2 


0 




This wing-flutter theory was applied to propellers by integrating 
the air forces over the blade as follows: For blades with constant 

chord, the velocity and hence l/k varies directly with radius; there- 
fore, the air forces must be integrated with respect to t]' which 
H + X 

equals . Since the elements of the determinant are obtained by 

H + L 

integrating with respect to t), it is advisable to set down the air-force 
terms in such a form that they are slLso functions of t) instead of tj ' : 

(5a) 


(5b) 


(5c) 


Therefore , 


q' = 


H + X 
H + L 


L 


q' = 


H + nL 
H + L 


At this point it appesirs most convenient to set up the integrals 
involving the air-force terms in the form of sumnations for use in a 
solution hy strip- analysis. 
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m + 1 


( 6 ) 


the blade Is equal to 


By taking the reference section at the tip, l/k at any point along 

^ h , ‘ 

— i — n ' — £ which would correspond to the resultant 
ktip" b , 

velocity for the condition of zero forwaiSi velocity. In forward flight 
the resultant velocity along the blade woulc^not vary linearly with radius 
and would be a function also of the forward velocity. For zero forward 
velocity, l/k according to equation (5c) becomes 


H + 


^tip S ^ 


h 

b 


iHL ‘ir 


(7) 


Equation (6) can then be written as follows: 




-2G + 2ip) ^ ^ 

/ b H + L ktip 


At] 


It should be noted that the aerodynamic coefficients F and G are 
related to the local values of l/k and thus vary along the blade radius 
also. For the purpose of strip analysisy ti must be measured to the 
center of each strip. 

Continuing the same proced\ire for each of the four determinant 
elements, equations (3), results in the following equations. Only the 
parts containing the air-force terms are sho^m as summations, because 
the mass terms can be Integrated mathematically for untapered blades. 




MCA TN 3357 


19 




20 


MCA TK 3557 


Equations (8) are substituted into equation (4) and solved -with the ■ 
final resTilt in the form of V/bc%, as; a function of The effect 

of centrifugal force on the static bending frequency can be computed as 
shown in reference 12. Centrifugal fdj:‘ce~also affects the static torsion 
frequency, and, for the present case, the_ same relationship as that used 
to correct the bending frequency was used as a first approximation for 
the corrected torsion frequency. The corrected ratio can be 

computed and plotted on the same graph with the flutter calculation. The 
intersection of the two curves yields Jthe theoretical flutter-speed 
coefficient for the given propeller. , _ 

The indicated theoretical flutter-speed coefficient is based on the 
torsional frequency corrected for cent^rlfugal force. In order to compare 
theory with experiment, the theoretical flutter-speed coefficient should 
be raised by the ratio of the corrected torsional frequency to the static 
torsional frequency. — 

Mode shapes of uniform untapered beams are presented in reference 12 
A method of obtaining mode shapes for nonunifoim beams and beams with 
concentrated masses is presented in ‘reference 13. 
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TABLE I 


tBSCmPPEOT <S FECESLLER-EIDHXKB IKDaLS 


[a* radlna ■ 1.21? 


Model 

3 iCA 

airfoil 

sactloa 

Bladfi 

■oterlal 

■TiS 

I. 

(ft) 

Ceps) 

f^qi 

(epa) 

CG 

(pare ant 
chord) 

EA 

(pereaat 

chco:d) 

1 /k 

{*) 


8 h 

Sd 

OJ 

(lb-ft2) 

BesasxfkB 

la. 

lb 

ifi-OflS 

1£-003 

allo^ 

Alntdnnn allqy 

0.167 

.167 

1.788 

1.370 

7.6 

12.6 

91.t 

122.0 

W .5 

W .5 

*5 

*5 

63 

63 

0.23 

.23 

0.011 

.011 

0.00* 

.00* 

36 


Ic 

lS-003 

Alunlima allc7 

.167 

.870 

32.1 

189.0 

W .5 

*5 

83 

.23 

.Oil 

.00* 

36 

17° twlBt, 

2 

i£.«q 3 

Alunlmii allCT 

.167 

1.788 

7.0 

82.5 

W .5 

50 

63 

.23 

.011 

.00* 

36 

3 a 

lfi-003 

Stael 

.167 

1.788 

7.2 

90.5 

* 8.5 

*7 

177 

.23 

.00* 

.001 

101 

root tc tip 


16-OC3 

Btaal 

.167 

1.788 

7.6 

95.0 

*8.5 

*6 

177 

.23 

.00* 

.001 

101 


3 c 

i£-0 C3 

Stael 

.167 

1.370 

12.8 

12 l ^.5 

*8.5 

*6 

177 

.23 

.00* 

.001 

Id 


3 d 

16-OC8 

Steel 

.167 

.870 

31.7 

202.3 

*8.5 

*6 

177 

.23 

.00* 

.OOL 

Id 



16-006 

Mepla 

.167 

1.788 

13 .i 

81 .2 

* 8.5 

*7 


.23 

.02 

.03 

12 


3 

15.:oo6 

A'Vnii'fwnj" a£[jo7 

'.250 

T.T88- 

•20:6‘ 

-176.0- 

-* 8.5 

** 

127 - 

-. 23 - 

-.OUr- 

,00t 

1*00 


6 

**16-003 

Alanlmm alloy 

.199 

1.750 

^3.7 

260.? 

*8.5 


76 

.23 

.011 

.00* 

---- 

Taper ratio ■ 0,57 

7 

16-006 

Maple with 

■ 

1.7^ 

11.3 

63.5 

37.5 

*2 

*6 

.32 

.02 

.03 

18 


8 

16-006 

hraae Inaart 
at 1aad^Tlg 

edge 

Ifc^la with 

.167 

1.788 

10.5 

56.2 

3*.0 

*0 

56 

.36 

.02 

.03 

20 


9 

16-006 

hraas tnaart 
a(t leading 
edge 
Maple 

.167 

1.788 

13.9 

73.0 

*8.5 

65 

28 

.23 

.02 

.03 

12 

Swept baiek 10® 

10 

16-006 

Heple 

,167 

1.7B8 

13.7 

78.6 

*8.5 

65 

28 

.23 

.02 

.C8 

12 

Swept back 20 ° 

U 

l£-oo6 

Maple 

.167 

1.708 

13.9 

96.9 

* 8.5 

62 

28 

.23 

.02 

.eg 

12 

Curred sweep- 

12 

16-006 

Kaple 

.200 

i.'Tsa 

30.0 

130.3 

* 8.5 

50 

26 

.23 

.02 

.03 


hack, 10° at tip 
Taper ratio - O.50 

13 

16-009 

Maple 

.167 

1.788 

19.0 

113.0 

* 8.5 

** 

*2 • 

.23 

.02 

.03 

*0 



•I/k cc«viit«d at p ■ O.OOC38 Bing par aibJje foot. 

-**Tapei«(l In tMcknaaa, t/Sb at root > 0,06 and t/9> at 0.9L •• 0.03. 
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TABLE U 

SCOPE OF PROPELLER-FLUTTER INVESTIGATIONS 


Parameters 

studied 

Range of 
values 

Models used 
for studies 

Figures in which 
parameters are 
evaluated 

Torsional stiffness 

12 to 101 lh-ft2 

la,3a,3h,4 

4 

Taper ratio 

0.50 to 1.0 

4,5,6,12 

5 

Blade twist (at tip) 

0° and 17 ° 

la, 2 

6 

! Length- chord ratio 

2.6 to 5-^ 

la,lh,lc,3a, 

3t,3c,3d 

7 

Density of operating 
medium 

0.0006 to 0.0021}- 
sl\ig/cu ft 

la 

8 

Section thickness ratio 

3 to 9 percent 
chord 

Is, 3a, 313,4, 13 

9 

Sweephack 

0° to 20° 

4,9,10,11 

10 

Section center-of- 
gravity location 

34.0 to 
48.5 percent 
chord 

4,7,8 

11 

Mach mzmher 

0 to 1.3 

5,6 

12 and 13 

Blade angle at 0.8L 

5° to 35° 

All 

3 to 13 








TABIE III 

CCMPABISON OF CLASSICAL-FnJTTER THEORT WITH EXPERIMEMTAL Dm 





y 

f\t\ 

ff=.\ 


Model 

Theory of 
reference 5 

(yA^^‘h)t±p 

rjM ju 1 j_ 1 rr\j UXUCUU. JT 

of reference 10 
(apTCndix) 

1 »t i iiimin_i_t_ii_MT._i_ 

wing theory of 
reference 7, 

vMt 

Theoretical 
representative section 
Cpercent blade length) 

Experimental data 

Classical 

flutter 

(v/baOtip 

Stall 

flatter 

C^A’^bcJtip 

2 ’ 

— 

4745 

■ 3.70 

75 

■■ 5.H ■ 

■ l.l4 

3a 

— 

7.16 

6.01 

76.5 

— 

1.24 


3-34 

3.a5 

2.68 

75 

3.4i 

l.l4 
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Figure 3*- Relation between lift coefficient and "blade-angle setting. 
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Figure 4.- Effect of torsional stiffness on flutter- speed 
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Figure 5*- Effect of tlade taper on flutter- speed coefficient 
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Figure 6.- Effect of "blade twist on flutter-speed coefficient. 
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(a) Models la, lb, and Ic. 

Figure J.- Effect of length-chord ratio on flutter-speed coefficient. 
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(b) Models 3a# 3b, 3c, and 3i. 
Figure 7.- Concluded. 
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Figure 9-- Effect of section thickness ratio on flutter-speed coefficient 
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Figure 10.- Effect of sweepback on flutter-speed coefficient 
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(a) Model 5, Po.8L = 25°. 


Figure 12.- Effect of Macli number on flutter- speed coefficient 
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(•b) Model 6, = 20° 

Figure 12.- Concluded. 
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Figure 13.- Effect of Mach, number on propeller flutter; model 6. 





